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Explicit current-dependent expressions for anisotropic longitudinal and transverse nonlinear mag- 
netoresistivities are represented and analyzed on the basis of a Fokker-Planck approach for two- 
dimensional single-vortex dynamics in a washboard pinning potential in the presence of point-like 
disorder. Graphical analysis of the resistive responses is presented both in the current-angle coordi- 
nates and in the rotating current scheme. The model describes nonlinear anisotropy effects caused 
by the competition of point-like (isotropic) and anisotropic pinning. Nonlinear guiding effects are 
discussed and the critical current anisotropy is analyzed. Gradually increasing the magnitude of 
isotropic pinning force this theory predicts a gradual decrease of the anisotropy of the magnetore- 
sistivities. The physics of transition from the new scaling relations for anisotropic Hall resistance 
in the absence of point-like pins to the well-known scaling relations for the point-like disorder is 
elucidated. This is discussed in terms of a gradual isotropizaton of the guided vortex motion, which 
is responsible for the existence in a washboard pinning potential of new (with respect to magnetic 
field reversal) Hall voltages. It is shown that whereas the Hall conductivity is not changed by pin- 
ning, the Hall resistivity can change its sign in some current-angle range due to presence of the 
competition between i- and a-pins. 
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I. INTRODUCTION 

The importance of flux-line pinning in preserving the 
superconductivity in a magnetic field has been generally 
recognized since the discovery of type-II superconductiv- 
ity. But till now the mechanism of flux-line pinning and 
creep in superconductors (and particularly in the high- 
T c superconductors (HTSC's)) is still a matter of contro- 
versy and great current interest, especially in the cases 
of strong competition between different types of pins. 

One of the open issues in the field is the influence 
of isotropic point-like disorder on the vortex dynamics 
in the anisotropic washboard planar pinning potential 
(PPP) for the case of arbitrary orientation of the trans- 
port current with respect to the PPP "channels" where 
the guiding of vortices can be realized. The importance 
of this issue may be substantiated by ubiquitous presence 
of point-like pins in those high- and low-T c superconduc- 
tors which were used so far for resistive measurements of 
the guided vortex motion 1-9 . 

The first attempt to discuss the influence of isotropic 
point-like disorder on the guiding of vortices was made 
by Niessen and Weijsenfeld 1 still in 1969. They stud- 
ied guided motion in the flux flow regime by measur- 
ing transverse voltages of cold-rolled sheets of a Nb-Ta 
alloy for different magnetic fields H, transport current 
densities J, temperatures T, and different angles a be- 
tween the rolling and current direction. The (H,J,T,a)- 
dependences of the cotangent of the angle (3 between the 



average vortex velocity (v) and the vector J direction 
were presented. For the discussion, a simple theoretical 
model was suggested, based on the assumption that vor- 
tex pinning and guiding can be described in terms of an 
isotropic pinning force F l p plus a pinning force F® with 
a fixed direction which was perpendicular to the rolling 
direction. The experimentally observed dependence of 
the transverse and longitudinal voltages on the magnetic 
field in the flux flow regime as a function of the angle a 
was in agreement with this model. 

Unfortunately, in spite of the correct description of 
a geometry of the motive forces of a problem (see 
below Fig. 1) it was impossible within the flux flow 
approach 1 to calculate theoretically the nonlinear (J, T, 
a)-dependences of the average pinning forces (F 1 ) 
and (Fp) which determine the experimentally observed 
cot/3(J, T, a) dependences. 

The nonlinear guiding problem was exactly solved at 
first only for the washboard PPP (i.e. for F l p = 0) 
within the framework of the two-dimensional single- 
vortex stochastic model of anisotropic pinning based on 
the Fokker-Planck equation with a concrete form of the 
pinning potential 10,11 . Two main reasons stimulated 
these theoretical studies. First, in some HTCS's twins 
can easily be formed during the crystal growth 2 " 5,8 . Sec- 
ond, in layered HTCS's the system of interlayers between 
parallel a6-planes can be considered as a set of unidirec- 
tional planar defects which provoke the intrinsic pinning 
of vortices 12 . 

Rather simple formulas were derived 11 for the ex- 
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perimentally observable nonlinear even(+) and odd(— ) 
(with respect to the magnetic field reversal) longitudi- 
nal and transverse magnetoresistivities pt ± (j,6,a,s) as 
functions of the dimcnsionless transport current density 
j, dimcnsionless temperature 9, and relative volume frac- 
tion < e < 1 occupied by the parallel twin planes di- 
rected at an angle a with respect to the current direction. 
The pjj^-formulas were presented as linear combinations 
of the even and odd parts of the function v(j, 9, a, e) 
which can be considered as the probability of overcoming 
the potential barrier of the twins 11 ; this made it possi- 
ble to give a simple physical treatment of the nonlinear 
regimes of vortex motion (see below item II. C). 

Besides the appearance of a relatively large even trans- 
verse resistivity, generated by the guiding of vortices 
along the channels of the washboard PPP, explicit ex- 
pressions for two new nonlinear anisotropic Hall resistiv- 
ities P|7 and were derived and analyzed. The physical 
origin of these odd contributions caused by the subtle in- 
terplay between even effect of vortex guiding and the odd 
Hall effect. Both new resistivities were going to zero in 
the linear regimes of the vortex motion (i.e. in the ther- 
moactivated flux flow (TAFF) and the ohmic flux flow 
(FF) regimes) and had a bump-like current or tempera- 
ture dependence in the vicinity of highly nonlinear resis- 
tive transition from the TAFF to the FF. As the new odd 
resistivities arose due to the Hall effect, their character- 
istic scale was proportional to the small Hall constant as 
for ordinary odd Hall effect investigated earlier 10 . It was 
shown 11 that appearance of these new odd Pi7j_ contribu- 
tions leads to the new specific angle-dependent "scaling" 
relations for the PPP which demonstrate the so-called 
anomalous Hall behavior in the type-II superconductors. 

Here we should to emphasize that the anomalous 
behavior of the Hall effect in many high-temperature 
and in some conventional superconductors in the mixed 
state remains one of the challenging issues in the vortex 
dynamics 5 ' 12,16 . The problem at issues includes several 
remarkable experimental facts: a) the Hall effect sign 
reversal in the vortex state with respect to the normal 
state at temperatures near T c and for moderate magnetic 
fields; b) the Hall resistivity "scaling" relation p± ~ p?, 
exists with 1 < f3 < 2, where p± is the Hall resistivity 
and /o | j is the longitudinal resistivity; c) the influence of 
pinning on the "Hall anomaly" and scaling relation. As- 
suming that the "bare" Hall coefficient an is constant, 
two different scaling laws have been derived earlier the- 
oretically for different pinning potentials 11,17 . Vinokur 
et al. have shown 17 that a scaling law p± = <5p 2 (where 

S = nanc 2 / B&o is the Hall conductivity, n = ±1, c 
is the speed of light, B is the magnetic field and $ is 
the magnetic flux quantum) is the general feature of any 
isotropic vortex dynamics with an average pinning force 
directed along the average vortex velocity vector. Later 
it was shown 11 that for purely anisotropic a-pins that 
create a washboard planar pinning potential, the form 
of corresponding "scaling" relation is highly anisotropic 



due to the reason that pinning force for a-pins is directed 
perpendicular to the pinning planes. If a is the angle 
between parallel pinning planes and direction of the cur- 
rent density vector j, then for a = the scaling law has 
the form p± = —n(otH /rj)p\\ (n is the vortex viscosity) 
which was interpreted previously 11 as a scaling law with 
0=1, whereas for a = n/2 the scaling relation is more 
complex 11 . The pj, as it is shown in this paper, can be 
presented as a sum of the three contributions with the 
different signs. The graphical analysis in Sec. Ill of this 
paper represents a some range of the (a, j)-values where 
the theory predicts a nonlinear change of the p^ sign. 

Let us consider another specific feature of the purely 
anisotropic guiding model 10 ' 11 . From the mathematical 
viewpoint, the nonlinear anisotropic problem, as solved 
in Ref . 1 1 , reduces to the Fokker-Planck equation of the 
one-dimensional vortex dynamics 13 because the vortex 
motion is unpinned in the direction which is parallel to 
the PPP channels. As a consequence, a critical current 
j c exists only for the direction which is strictly perpen- 
dicular to the PPP channels (a = 0); j c (a) = for any 
other direction (0 < a < ir/2). However, the measure- 
ments of the magnetoresistivity show 1-8 that j c (oi) > 
for all a (although j c (a) may be anisotropic). So, in 
spite of some merits of a model with a washboard PPP, 
which was the first exactly solvable stochastic nonlinear 
model of anisotropic pinning, it cannot describe the j c - 
anisotropy of the experimentally measured samples. 

Due to this reason later it was suggested 14 ' 20 another 
simple model, which demonstrates this j c -anisotropy for 
all a on the basis of the bianisotropic pinning potential 
formed by the sum of two washboard PPP's in two mutu- 
ally perpendicular directions. In contradistinction to the 
nonlinear model with uniaxial PPP 11 , this bianisotropic 
nonlinear model predicts a j c (a)-anisotropy and relates 
it to the guiding anisotropy, describing the appearance 
of two step-like and two bump-like singularities in the 
and p^ ± (Hall) resistive responses, respectively. Al- 
though several proposals to realize experimentally this 
bianisotropic model were discussed so far 14 , the corre- 
sponding experiments, however, are still absent. 

At the same time, the experimental study of vortex dy- 
namics in the PPP is always accompanied with a presence 
of a certain level of point-like disorder. So, as far as the 
analysis of existing experimental data is concerned, none 
of the present theoretical studies in the limiting cases of 
purely anisotropic or isotropic pinning arc sufficient. The 
more general approach is needed. 

The objective of this paper is to present results of a 
theory for the calculation of the nonlinear magnetoresis- 
tivity tensor at arbitrary value of competition between 
point-like and anisotropic planar disorder for the case of 
in-plane geometry of experiment. This approach will give 
us the experimentally important theoretical model which 
demonstrates the j c -anisotropy for all a and predicts a 
nonlinear change of the pj sign at some set of parameters 
(without change of the Hall conductivity) due a compe- 
tition of the washboard PPP and a point-like disorder. 



3 



The organization of the article is as follows. In Sec. 
II we derive main results of the i + a pinning problem 
and consider two main limiting cases of purely a- or i- 
pinning. In Sec. Ill we represent the graphical analysis 
of different types of nonlinear responses, in particular, 
the (J, a) graphs of the pjj^ magnetoresistivities and the 
resistive response in a rotating current scheme. In Sec. 
IV we conclude with a general discussion of our results. 

II. MAIN RELATIONS 

A. Formulation of the problem. 

The Langevin equation for a vortex moving with ve- 
locity v in a magnetic field B = nB (B = |B|, n = nz, 
z is the unit vector in the z-direction and n = ±1) has 
the form 

rjov + na H v xz = F L + F p ] + F p + F th , (1) 

where = n(Q /c)j x z is the Lorentz force ($ 1S the 
magnetic flux quantum, c is the speed of light), F° = 
— VU p (x) is the anisotropic pinning force (U p (x) is the 
washboard planar pinning potential), F^ is the isotropic 
pinning force, induced by uncorrelated point-like disorder 
, F t h is the thermal fluctuation force, 770 is the vortex 
viscosity, and an is the Hall constant. 

For purely isotropic pinning (i.e. for F p = 0) Eq. ^ 
was earlier solved 17 for F t h = 0, using the fact that 

f; = -m(v)v, (2) 

where rji(v) is velocity-dependent viscosity and v = |v|. 

Below we will show (see Eq. (8) and item D of Sec. 
II), that the solution, obtained in Ref. 17, can be pre- 
sented in terms of the probability function of overcoming 
the effective current- and temperature-dependent poten- 
tial barrier of isotropic pinning Vi(Fi), which is simply 
related to Tji(v). 

In the absence of point-like disorder (i.e. for F l p = 
0) Eq. ([1]) was reduced to the Fokker-Planck equation, 
which was solved 10,11 , assuming that the fluctuational 
force F t h(t) is represented by a Gaussian white noise, 
whose stochastic properties are assigned by the relations 

(FthM) = 0, (FthAt)Fth,j{t')) = 2T m 6 t3 5(t - t'), (3) 

where T is the temperature in energy units. 

In what follows we derive the solution of Eq. ([T]) , using 
for F p the assumption ([2]), which reduces Eq. (fTJ) to the 
equation 

!)v + M H vxz = F i + F;+r i/l , (4) 

where r\ = rj{v) = r/o + Wi(v). Using the result of Ref. 11, 
the selfconsistent solution of the Eq. ([4} can be repre- 
sented as 

r,{v){v x )=F aVa {F a )/{\ + l 2 ), 

(5) 

V(v)(v y ) = F Ly + neF a v a {F a )/{l + i 2 ), 




FIG. 1: System of coordinates xy (with the unit vectors x 
and y) associated with the PPP planes and the system of 
coordinates x'y' associated with the direction of the current 
density vector j; a is the angle between the channels of the 
PPP and j , P is the angle between the average velocity vector 
of the vortices v and the vector j; Fl is the Lorentz force; 
< Fp > and < Fp > are the average isotropic and anisotropic 
pinning forces, respectively, Fi is the average effective motive 
force for a vortex. Here for simplicity we assume e = 0. 

where v a {F a ) is the probability of overcoming the PPP 
under the influence the effective moving force F a = 
Flx — nlF^y, Fl x and F^ y are the Lorentz force com- 
ponents acting along the vector x and y, respectively, 
I = eZ{v), e = an /r]o and Z(v) = i]o/r](v) with an obvi- 
ous condition < Z[v) < 1. Eqs. ([5]) can be rewritten 
as 

ri(v)(v)=F x , (6) 

where Fj x and Fj v are corresponding right-hand parts of 
Eqs. ([H]). From Eq. ^ we have 

r ] {v)v = F I , (7) 

where Fj = (Ff x + Ff y Y^ 2 and we omitted for simplicity 
the symbol of averaging for v. Then from Eq. ([7]) follows 
that v — v(Fj) and thus it is possible to represent f]i(v) 
and Z(v) in terms of Fj\ Tji(v) = r/i[v — v(Fi)] = fji(Fi) 
and 

Z(v) = Z[v = v(F I )] = u i (F I ). (8) 

Here Vi{Fj) has a physical meaning of the probability 
to overcome the effective potential barrier of isotropic 
pinning under the influence of effective (u-dependent 
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through the e-dependence) force Fj. Then in terms of 
the Vi{Fj) Eq. ^ takes the selfconsistent form 

% v = ViiF^Fz, (9) 

which can be highly simplified for a small dimension- 
less Hall constant (e <C 1). Really, in this limit 
e = eVi(Fi), where Fj = Fj(e = 0), and the right- 
hand part of the Eq. ([6]) becomes v- independent, i.e. is 
represented only in terms of the known quantities. Just 
in this limit all subsequent results of the paper will be 
discussed. 



B. The nonlinear resistivity and conductivity 
tensors 

The average electric field induced by the moving vortex 
system is given by 

E = (l/c)B x v = n(B/c)(-v y x + v x y), (10) 

where x and y are the unit vectors in x- and y-direction, 
respectively. 

From formulas ([9]) and (flQ)) we find the dimension- 
less magnetoresistivity tensor p (having components mea- 
sured in units of the flux-flow resistivity p/ = QoB/r/oc 2 ) 
for the nonlinear law E — p(j)j 

p — ( P xx P x y J 

V Pyx Pyy J 

(11) 

/ Vi (Fi) -nevf(Fi)p a (F Lx ) \ 
\ nev?(Fi)v a {F Lx ) i/<(*»i/„(.F„) J ' 

The conductivity tensor a (the components of which 
are measured in units of 1/p/), which is the inverse of 
the tensor p, has the form 

g. — I a xx a xy \ _ 
\ a yx a yy J 

(12) 

/ [MFiT 1 rie \ 

\ -ne [vi{F I )v a (F a )Y 1 ) ' 

From Eqs. (fTTj) and (TT2]) we see that the off-diagonal 
components of the p and a tensors satisfy the Onsager 
relation (p xy = —p yx in the general nonlinear case and 
(j xy = —<J yx ). All the components of the p-tensor and the 
diagonal components of the cr-tensor are functions of the 
current density j through the external force value Fl, the 
temperature T, the angle a, and the dimensionless Hall 
parameter e. For the following (see item E.2 of Sec. II) 
it is important, however, to stress that the off-diagonal 
components of the a (i.e. the dimensional Hall conduc- 
tivity terms <5 = ne/ pf) are not influenced by a presence 
of the i- and a-pins 16 . 

The experimentally measurable resistive responses re- 
fer to a coordinate system tied to the current (see Fig. 1). 



The longitudinal and transverse (with respect to the cur- 
rent direction) components of the electric field, En and 
E±, are related to E x and E y by the simple expressions 

E\\ = E x sin a + E y cos a, 

E±_ = —E x cos a + E y sin a. 

Then according to Eqs. (TlU)) . the expressions for the 
experimentally observable longitudinal and transverse 
(with respect to the j-direction ) magnetoresistivities 

= E\\/j and p± = Ej_/j have the form: 

p\\ = Pxx sin 2 a + pyy cos 2 a, 

(14) 

P± = Pyx + (pyy - Pxx) SU1 a COS a. 

Note, however, that the magnitudes of the P\\,±, given 
by Eqs. (fl"4"l) , are, in general, depend on the direction of 
the external magnetic field B along z axis due to the ne- 
dependence of the Fj and F a forces in arguments of the i>\ 
and v a functions, respectively. In order to consider only 
n-independent magnitudes of the p\\- and p^ -resistivities 
we should introduce the even(+) and the odd(— ) with re- 
spect to magnetic field reversal (p^ = (p(n) ± p(—n))/2) 
longitudinal and transverse dimensional magnetoresistiv- 
ities, which in view of Eqs. (jT4j) have the form: 

pj = p/[sin 2 a + v a (F Lx ) cos 2 a\Vi(Fi), 

p^ = p s [[siv? a + v a {F Lx )cos 2 a]vr{F!)+ (15) 

v^F^v- (F a ) cos 2 a]. 



Px = -pfVi(Fi)[l - v a (F Lx )] sinacosa, 
Pl = p f [nev a (F Lx )v}(F l ) + 

(16) 

{v-{F a ) Vi {F^ - ur(F z )[l - v a (F Lx )}}x 

sin a cos a] . 

Here v~ are the odd (v~ = v~(n) = (v(n) — v(— n))/2) 
components of the functions Vi(Fj) and i> a (F a ), and for 
v~(F a ) we have the expansion in terms of e << 1: 

v- ~ -nevi(Fi)F L y[dv a (F Lx ) / dF Lx ]. (17) 

Eqs. (jX5j)-(jT6j) are accurate to the first order in e << 1 
and contain a lot of new physical information, which will 
be analyzed below (see item E of Sec. II). However, be- 
fore this analysis it is instructive to discuss in short the 
main physically important features of two main limit- 
ing cases of purely anisotropic a-pinning and isotropic 
i-pinning, which follow from Eqs. (115p - (116|) . when i/j = 1 
or v a = 1, respectively. 
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C. Anisotropic a-pinning. 

Setting i/j = 1 we obtain rather simple formulas, 
which were derived firstly 11 for the experimentally ob- 
servable nonlinear even and odd longitudinal and trans- 
verse anisotropic magnetoresistivities p^ j_(j,9,a,s a ) as 
functions of the transport current density j, dimension- 
less temperature 9 and relative volume fraction < e a < 
1, occupied by the parallel twin planes, directed at an 
angle a with respect to the current direction: 

P\\ a = Pf[v£ cos 2 a+sin 2 a], p+ a = p f (v+-l) sin a cos a, 

(18) 

P\\a = Pf v a cos 2 a, p1 a = p f [nev+ + v~ sin a cos a] . 

(19) 

Here v a = v a {F) is considered as the probability of 
overcoming the potential barrier of the washboard PPP 
in the x-direction under the influence of the effective 
force F = Fl x — neF^y 11 . This ^-function describes 
an essentially nonlinear transition from the linear low- 
temperature thermoactivated flux flow (TAFF) regime 
of vortex motion to the ohmic flux flow (FF) regime. It 
is a step-like function of j or 9 for a small fixed temper- 
ature or current density respectively (see Figs. 4, 5 in 
Rcf. 11). 

It follows from Eqs. (18)-(19) that for a ^ 0,tt/2 the 
observed resistive response contains not only the ordi- 
nary longitudinal p^ a (a) and transverse p~J_ a (a) magne- 
toresistivities, but also two new components induced by 
the pinning anisotropy: an even transverse p\_ a {ot) and 
an odd longitudinal component p7,{a). The physical ori- 
gin of the p±„(ct) (which is independent of e) is related 
in an obvious way with the guided vortex motion along 
the " channels" of the washboard pinning potential in the 
TAFF regime. On the other hand, the component p~ a {ct) 
is proportional to the odd component v~ , which is zero 
at e = and has a maximum in the region of the nonlin- 
ear transition from the TAFF to the FF regime at e ^ 
(see Figs. 6, 7 in Ref. 11) The (j, 6*)-dependence of the 
odd transverse (Hall) resistivity p± a (j,0) has contribu- 
tions both, from the even vj w v a and from the odd 
v~ components of the u a (j, 0)-function. Their relative 
magnitudes are determined by the angle a and the effec- 
tive Hall constant e. Note, that as the odd longitudinal 
p7 and odd transverse pJ Q magnetoresistivities arise by 
virtue of the Hall effect, their characteristic scale is pro- 
portional to e << 1 (see Eqs. (19)). 

The appearance of these new odd Hall contributions 
follows from emergence of a certain equivalence of xy- 
directions for the case, where a guiding of the vortex 
along the channels of the washboard anisotropic pinning 
potential is realized 18 at a ^ 0, 7r/2 and leads to the new 
specific angle-dependent "scaling" relations for the Hall 
conductivity 11 for the case etana << 1 

ne = (p1 a - p« tana) cos 2 a/ (p+ - p f sin 2 a). (20) 



Here the dimensionless Hall constant e << 1 is uniquely 
related to three experimentally observable nonlinear re- 
sistivities Piu) Piu) Pla) an d the "scaling" relation (|20|) 
depends on the angle a. This relation differs substan- 
tially from the power-law scaling relations, obtained in 
the isotropic case 17 (see below). In the particular case 
a = we regain the results 10 , specifically e = p± a / pt a 
(in Ref. 10 e = p±/ p\\), i.e. a linear relationship between 

Pla and Pta' 

Eq. (20) may be represented in another form 

Pla ( a ) =5v a {a)pf 2 -P\\ a {a) tana (21) 

which is more suitable for considering scaling relations 
in longitudinal (a = 7r/2) and transverse (a = 0) LT- 
geometries of experiment 11 . In these geometries second 
term in the right hand side of Eq. (21) is zero and we 
obtain that 

Pl a = Hpl la f, (22) 
5(a = ir/2)sS L = 6va(0,B)\ 

(23) 

S(a = 0) = S T = 5/u a {j,e). 

From Eqs. (22)- (23) follows that 5 may be interpreted 
as an effective Hall conductivity in LT-geometries which 
is suppressed for a = ir/2 (5l < 6) and enhanced for 
a = (St > 5) in comparison with a bare Hall con- 
ductivity S. The physical reason for this influence of v a - 
function on the S behavior in LT-geometries is simple. 
Namely, it appears as a result of the fact that in the case 
of anisotropic pinning the driving force F = FL x — neFL Vl 
which determines the probability of overcoming the po- 
tential barrier ( and therewith also determines the mag- 
nitude of the component of the vortex velocity perpen- 
dicular to the channels of the PPP), is the sum of two 
forces. The first of these is the transverse component 
of the Lorentz force, Fl x = i^cosa, and the other is 
the transverse component of the Hall force Fh — neFLy 
which is proportional to the longitudinal (relative to the 
PPP planes) component of the velocity of guided vortex 
motion. This second force Fjj, which changes its sign 
(relative to the sign of Fl ) upon reversal of the sign of 
the external magnetic field, is the reason for appearance 
of new, Hall-like in their origin, i/~-terms in the formulas 
for the resistive responses in Eqs. (19). 

Returning to the physics of suppression and enhance- 
ment of the S in LT-geometries we should keep in mind 
that only longitudinal component of the vortex velocity 
(with respect to the current direction) vi is responsible 
for the appearance of the transverse Hall voltage. Thus, 
in L-geometry vi and S are suppressed by PPP-barriers, 
whereas in T-geometry vi is not influenced by them and 
5 looks like enhanced quantity. On the contrary, the be- 
havior of the transverse component of the vortex veloc- 
ity vt, which determines the longitudinal voltage, in LT- 
geometries is opposite. 



6 



D. Isotropic i-pinning. 

For this case we put v a = 1 and from Eqs. (fl"5"]) - (fl"6]) 
follows that 

P\\ = P\\i = PfVi(F L ), p±= p±i= pfnev?(F L ), (24) 



where F L = Fi{v a = 1) =| F L |. From Eqs. ([24]) the 
well-known scaling relation piu ~ (p±i) 2 , derived firstly 
in Ref. 17, follows. Note that = p7. — in this case, 
i.e. nonlinear resistive response is isotropic. 

E. Competition between a- and i-pinning. 



Equations ([15 ]) - ([16]) for the magnetoresistivities p« ± 
at arbitrary value of competition between point-like and 
anisotropic planar disorder for the in-plane geometry of 
experiment can be represented in a more suitable form, 
if we take into account Eqs. ([18])-([19]) and ([24]) : 

p+ = v l (F i )-p+ p+ = v l (F l )-p+ a , (25) 



P\\ = v i P \ 



(26) 



p ± = p f nev a vf + Pf{v a v; t -v i [1 - v a }} sin2a/2. (27) 

Here Vi(Fi) is the probability function vi of anisotropic 
argument F 2 = [Fl x vl(F Lx ) + F^} 1 / 2 , the magne- 
toresistivity pt± a an d the v a = ^ Q (i^L:r)-furictions in 
Eqs. ([25]) - ([27]) are the same as those in item C of Sec. 
II; vr = urlF^n)] and F T (n) = [F% + F 2 Lx v 2 a (F a ) + 
2neVi{Fi)FL X FL y v a {\ — i^)] 1 / 2 . It is easy to check, that 
previous results of items C and D of Sec. II follow from 
Eqs. (|2l))) - J27p in the limits of purely anisotropic (i.e. for 
Vj, = 1, za~ = 0) and isotropic (i.e. for v a = 1, v~ = 0) 
pins. 

In this subsection it must be suffice to discuss in short 
the main physically important features of these equa- 
tions. First of all, the magnetoresistivities pi ± can be 
found, if the v a - and v^- functions are known. Moreover, 
the converse statement is also valid: it is possible to re- 
construct these functions from (j, 8, independent resis- 
tive measurements, using only Eqs. (|25|) . where the Hall 
terms are ignored. Eqs. ([2"6"]) and (|27p. which arise due 
to the Hall effect, have a rather complicated structure, 
which reflects a more pronounced competition between 
isotropic and anisotropic disorder in the Hall-mediated 
resistive responses. Let us outline the main new physical 
results, following from Eqs. ([25]) - ([27]) . 



1. Point-like disorder and vortex guiding. 

For the discussion of the influence of point-like pins 
on the guiding of vortices in the anisotropic pinning po- 
tential it is sufficient to analyze Eqs. E5]) , Whereas for 



the purely anisotropic pinning (i/j = 1) a critical cur- 
rent density j c exists only for direction, which is strictly 
perpendicular to the PPP (a — 0) and j c (a) — for 
any other direction (0 < a < tt/2) due to the guiding 
of vortices along the channels of a washboard potential, 
in Eqs. (|25|) the factor Vi(Fi) ensures that an anisotropic 
critical current density j c (a,9) exists for arbitrary an- 
gles a. 

It is interesting, however, to note, that the angular de- 
pendence of the ratio p±/p\\, which determines the angle 
(3 between j and v for a-pins in Ref. 11, according to the 
relation 



cot/3 = -^a = 



1 - V a 



P|| a tan a + v a cot a 



(28) 



is not influenced by the isotropic disorder, because factor 
Vi(Fi) in Eqs. ([25f vanishes from Eq. (j28|) . Physically it 
means, that character of anisotropy in the case of com- 
petition between i- and a-pinning is determined only by 

< F; >= [v a (F Lx ,0) - l]F Lx , (see Fig. 1), i.e. by the 
average pinning force of the PPP. Isotropic pins influ- 
ence only the magnitude of the average v-vector, because 

< F* >|| v || Fj. So, the polar resistivity diagram p(a), 
which can be measured experimentally 5 , is influenced by 
point-like pins, because from Eqs. (fTTj) follows, that 



p(a) = p f [p\ x sin 2 a + pi cos 2 a] 



1/2 



PfVi(Fi) (sin 2 a + v 2 cos 2 a) 1 / 2 . 



(29) 



2. New Hall voltages and scaling relations. 

As it follows from Eqs. ([26]) -([27]), the odd longitudinal 
p7 and transverse magnetoresistivities contain terms 

with the ^"-function. They possess a highly anisotropic 
current- and temperature-dependent bump-like behavior. 
They tend to zero in the linear regime of vortex motion. 
For a = 0, 7r/2 these new terms disappear, because vf = 
v~ = at these limits. As it was in the case of purely 
a-pinning (see item C of Sec. II), the appearance of these 
new odd Hall contributions follows from the emergence 
of a certain equivalence of xy-directions due to a guiding 
of vortices along the channels of the washboard pinning 
potential for the case with a ^ 0, 7r/2. Note also, that p7 

includes two terms with similar signs, whereas in there 
are terms with opposite signs. The latter can give rise to 
the well-known sign change in the (J, 8, H )-dependence 
of the Hall resistivity below T c 12 . 

From Eqs. ([25]) - ([27]) new anisotropic "scaling" rela- 
tions for the dimensionless Hall constant e can be de- 
rived. For this purpose we exclude z/~ from Eqs. ([23]) - 
([2"T]) . for v~ use Eq. ([T7]) ; and after some algebra in the 
limit e • tana << 1 we have: 

2p]| ■ pt + pf sin 2a(l - v a )vip7 

ne = r — y. (30) 

[2v a pT - sin 2a • PfV z F Ly v a \pfvf 
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It is easy to check that from Eq. (|30| follows scaling 
relations S — ne/pj = p±i/(p\\i) 2 (for i-pins at v a = 1) 
and Eq. (|20|) (for a-pins at Vi = 1) 

As it follows from Eqs. (25) and (27), just the same 
"scaling" relations as given by Eqs. (22) and (23) for a- 
pins, exist also for (i+a)-pins (with a replacement of cor- 
responding a-resistivities in Eq. (8) by pjp~J_ and pfp^). 
Physically it follows from the fact that point-like disor- 
der does not change the angular dependence of the ratio 
P~±al P\\ a > wn i cn determines the angle (3 between j and av- 
erage velocity vector (v) for a-pins 11 , and influences only 
the magnitude of (v) 16 . 



III. GRAFICAL ANALYSIS OF NONLINEAR 
REGIMES. 

A. Pinning potential and ^-function behavior. 

In order to analyze different types of nonlinear 
anisotropic (j, 9, a)-dependent magnetoresistivity re- 
sponses, given by formulas ([25|) -([27 | . we should bear 
in mind that these responses, as is seen from formula 
(fTTj) . are completely determined by the (j, 6*)-behavior 
of the functions v a (F a ) and ^(Fj), having a sense of 
the probabilities to overcome the effective potential bar- 
riers of the a- and i-pins, respectively. A simple ana- 
lytical model for the calculation of the (j, #)-dependent 
^-functions was given earlier 11,13,20 . We will use for both 
Vi and v a functions the one-dimensional periodic pinning 
potential U p (x) (see Fig. 2), which has a simple analytical 
form 11,20 : 



Up * 





F p x, 



< x < b, 



U p (x) = { F p {x - 26), bs^x^ 2b, (31) 



0. 



2b < x ^ h, 



where F p is the pinning force (F p — Uo/b, where Uo > 
is the depth of the potential well and 2b is the width 
of the well). This form of U p (x) allows to define as the 
properties of a given pinning center (by the parameters 
Uo and b) , as well as the density of such centers (by the 
parameter e = 2b/ h, where h is the period of the U p (x)). 

Calculation of the v(j, 9) function on the basis of the 
pinning potential, given by Eq. (|3ip . was done 11 and can 
be represented here in the form 11 



2/(/ 2 - I) 2 



2/(/ 2 -l)(/ 2 -l + e)-e6 1 G' 
where G = {(3/ 2 + 1) cosh(//0£) + (f 2 - 1) cosh[(/(l 
2e))/(9e)} - 2/(/ - 1) cosh[/(l - e)/9e - (1/6)} - 
2f(f + 1) cosh[(/(l - e)/9e + 1/9}}/ sinh(//fe). 



(32) 



■U<> 



2b 





FIG. 2: Model pinning potential U p (x): h is the period of the 
potential, 26 is the width of the potential well, Uo is the depth 
of the potential well, e = 2b /h characterizes the concentration 
of the pinning planes. 



Here and below we have for the time being dropped 
the indices a and i from the physical quantities pertain- 
ing to pinning potentials U pa and U P i and formula (|3"2"]) 
describes equally the pinning on both potentials. For 
convenience of qualitative analysis of the formulas fol- 
lowing dimensionless parameters were used: / = Fb/Uo 
is the effective motive force, which specifies its ratio to 
the pinning force F p — Uo/b, 9 = T/Uq is the tempera- 
ture. 

The effect of the external force F acting on the vortices 
consists in a lowering of the potential barrier for vortices 
localized at pinning centers and, hence, an increase in 
their probability of escape from them. Increasing the 
temperature also leads to an increase in the probability 
to escape of the vortices from the pinning centers through 
an increase in the energy of thermal fluctuations of the 
vortices. Thus the pinning potential of a pinning center, 
which for F, T — > leads to localization of the vortices, 
can be suppressed by both an external force and by tem- 
perature. 

A detailed quantitative and qualitative analysis of the 
behavior of v(f, 9, e) as a function of all the parameters 
and its asymptotic behavior as a function of each are 
described 11 . Here we will pay particular attention only 
to the typical curves of v{f, 9, e) as a function of the pa- 
rameters / and 9, which describe the nonlinear dynamics 
of the vortex system as a function of the external force 
acting on the vortices in the direction perpendicular to 
the pinning centers and as a function of temperature (see 
Figs. 4 and 5 in Ref. 11). As we see from those figures, 
the form of the v(f) and v(9) curves is determined by 
the values of the fixed parameters 9 and /. The mono- 
tonically increasing function v(f) reflects the nonlinear 
transition of the vortex motion from the TAFF to the FF 
regime with the increasing external force at low temper- 
atures (T <C Uo), while at high temperatures (T 3> Uq) 
the FF regime is realized in the entire range of variation 
of the external force (even at small forces) because of the 
effect of thermal fluctuation on the vortices. The mono- 
tonically increasing function v(9) reflects the nonlinear 
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transition from a dynamical state corresponding to the 
value of the external force at zero temperature to the FF 
saturation regime. The width of the transition from the 
TAFF to the FF regime on the v(f) and v(9) curves de- 
pends on substantially different on the increasing of the 
parameters 9 and /, respectively. Namely, with increas- 
ing 9 the function v(f) shifts leftward and becomes less 
steep (see Fig. 4 in Ref. 11). That is, the higher the tem- 
perature, the smoother the transition from the TAFF to 
the FF regime and the lower the values of the external 
force, at which it occurs. With increasing / the v(9) 
curve also shifts leftward, it becomes steeper (see Fig. 5 
in Ref. 11). Consequently, the greater the suppression of 
the potential barrier of the pinning center by the exter- 
nal force, the sharper the transition from the TAFF to 
the FF regime and the lower the temperature at which 
it occurs. 

These graphs will be needed later on when we will dis- 
cuss the physical interpretation of the observed guiding- 
depended resistive responses. We also note that the de- 
pendence of the probability function v(e) on the concen- 
tration of pinning centers decreases monotonically from 
the value i/(Q) = 1, which corresponds to the absence 
of pinning centers, and that it becomes steeper with de- 
creasing fixed parameters / and 9, owing to the growth 
of the probability density for finding the vortices at the 
pinning centers with decreasing temperature and exter- 
nal force. 



B. Dimensionless form of the p„ ± -responses. 

Let us turn to the dimensionless parameters by which 
one can in general case take into account the difference 
of the potentials U a and Ui~ specifically, the difference 
of their periods h ai hi, the potential well depths Uo a , 
Uoi and the width b a , bi. We introduce some new pa- 
rameters: e = (SaEi) 1 ^ 2 is the average concentration 
of pinning centers, Uq = (IW^oi) 1 ^ 2 i s the average 
depth of potential well, k — (ei/ea) 1 ^ 2 = {habi/hiba) 1 / 2 , 
and p = (Uoa/Uoi) 1 / 2 , where the parameters k and 
p are measures of the corresponding anisotropies. The 
temperature will be characterized by new parameters: 
9 a = pT/Uo = T/U 0a and 9 t = (l/p)T/U = T/U 0i , 
which are the ratio of the energy of thermal fluctuations 
of the vortices to the average potential well depth U^ a 
and Uoi, respectively. 

The current density will be measured in units of j c = 
cUo/<&oh, where h = (hahi) 1 / 2 . Then the dimensionless 
parameters f a and /j, which specify the ratio of the exter- 
nal forces F a and F{ to the pinning forces F pa = Uo a /b a 
and F p i — U$i/bi (y a and Vi are the even functions of their 
arguments), we denote as f a = F a /F pa and = Fi/F pl . 
The values of the external force F, at which the heights 
of the potential barriers Uo a and Uoi vanish at T = 
correspond (at a = and a — tt/2) to the critical 
current densities j ca = qj c and jd = j c /q respectively, 
where q = p/n. In general case of nonzero temperature 



and < a < n/2 it is possible to consider the angle- 
dependent crossover current densities j C a{ot) an d j C i{a) 
(see below) which correspond to change in the vortex 
dynamics from the TAFF regime to a nonlinear regime. 
The condition, that determines the temperature region, 
in which the concept of critical current densities is phys- 
ically meaningful is ^ T <C Uq, because for T > Ua 
the transition from the TAFF to the nonlinear regime 
is smeared, and the concept of critical current loses its 
physical meaning. 

It is possible now to rewrite Eqs. (jT5])-(|T5|) in the di- 
mensionless form in order to represent them as functions 
of j, 9, a at given values of parameters e, e, q, k. 

p{ = ^(/ 4 )[sin 2 a + v a {f a ) cos 2 a)], (33) 



Pi = -n{fi)[l - va(fa)} sin 2a/2, (34) 



P\\ = v i {fi)[sm 2 a + v a (f a ) cos 2 a)] + v a (f a )vi(fi) cos 2 a, 

(35) 



p ± = nev a (f a )v 2 (f t ) + 

(36) 

+K"(/aM(/i) - - Va(fa)}} SUl 2a/2, 



where fa=jq 1 cosa, (37) 



fi= jq(sm 2 a + v 2 a cos 2 a) 1 ' 2 , (38) 

and /„ = jq [ei/i (fi ) sin a + n cos a] , 

ft = Msin 2 a + vl(fa) cos 2 a) - nev i {f i )v a {f a ) x 

x[l- i y Q (/ a )]sin2a} 1 /2. 

Here 

Va(fa) = V a {f a ,6 a ,e a /n), Va{fa) = V a {fa, a , E a / '«) , 

Vi{fo) = Pi(fi,6i,£iK), Vi{fi) = Vi(fi,6i,EiK), 

and vfifi) = {uilUn)] ± Ui[fi(-n)]}/2, 
vtifa) = Kfe] ± v a [U-n)]}/2. 

In Eqs. (33)-(38) we also denoted Vi(fi) = vf{fi) and 
Va(fa) = v£(f a ) for simplicity. 
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Before following graphical analysis of the p^ ± depen- 
dences given by Eqs. (|33|) -([36 |) . we should point out the 
magnitude of some parameters which will be used for pre- 
sentation of the graphs. It is important to remind here 
that the parameter q determines the value of anisotropy 
between Vi and v a critical current densities, whereas the 
parameter k describes the anisotropy magnitude of the 
width of nonlinear transition from the TAFF to the FF 
regime for vi and v a function. More definitely, if q > 1, 
then j ca = qj c > jd = j c /q and influence of the i-pins 
on the vortex dynamics decreases with q-increasing. For 
q < 1 the situation is opposite and anisotropy effects may 
be fully suppressed with ^-decreasing. So, for the obser- 
vation of pronounced competition between i- and a-pins 
q « 1 should be taken. 

The temperature dependences of the pjj~(a) at small 
current densities under conditions of the presence both 
isotropic and anisotropic pinning potential were stud- 
ied experimentally 8 . Arrhcnius analysis of these de- 
pendences within the frames of suggested here theo- 
retical approach have shown that for the samples 8 the 
U 0a = 4031K, U 0i = 1568K, b a = 400 nm, b t = 2000 
nm at T w 8K. Then for these samples q « 1.6, K « 0.5, 
w 0.003. It was also pointed out 8 that the best fitting of 
the experimental and theoretical curves was established 
for bijb a — 15, from which follows k ~ 0.25. So for all 
graphs below we used q — 1.6, k — 0.25, 9 = 0.003, 
e = 0.01 and if it is not pointed out specially, e a — 1 and 

Ei = 0.1. 

Note also that for the even longitudinal resistivity 
and the even transverse resistivity p\_ for a small 
Hall effect, terms proportional to e « 1 are absent 
(see Eqs. (l33l) - ((3"4")) ) and only contributions describing 
the competition between isotropic pinning and nonlin- 
ear guiding effect on the PPP in terms of the even j/j and 
u a functions are presented. 



C. Graphical analysis of current-angular 
dependences. 

1. (j, a) -presentation of v a and Vi. 

In order to discuss graphical (j, a)-behavior of the re- 
sistive responses we will use v a and Vi functions of their 
arguments f a and respectively, in the form given by 
Eq. (|3l?|) . Then these functions are, as a corresponding 
^-function 11 , the step-functions in j (at fixed 6) or in 
6 (at fixed j). For every of the i/-functions it is useful 
to determine the "crossover current densities" jci{ca) and 
Jca (a), as those which correspond to the middle point 
of a sharp step-like nonlinear transition from the TAFF 
to the FF regime. As it follows from Eqs. l[37] l -([55 |) . we 
can present /„ and as f a = j/j ca (a) with j ca (a) = 
q/cosa, and f+ = j/jd(a) with jd(a) rs l/qcosa for 
a <C 7r/4 and j C i(a) ~ l/v a qcosa for tan 2 a v%(j,a); 
jci{ot) w 1/gsina for a > 7r/4 because Eq. can 




FIG. 3: The current-angle dependence of the anisotropic 
probability function v a {j,ot). In all following graphs the pa- 
rameters q = 1.6, k — 0.25, 9 = 0.003, e = 0.01, e a = 1, and 
Si =0.1 (unless otherwise stated). 




FIG. 4: The current-angle dependence of the average effective 
motive force for a vortex fi(j, a). 

be presented in two equivalent forms, namely fi — 
jq cos ai/tan a 2 + = jqsina^/l + (f a /tana) 2 . 

The behavior of u a (j, a) function (see Fig. 3) is rather 
evident from the f a (j,ct) behavior. Namely, for all a ^ 
7r/2 (i.e. for f a ^ 0)) the v a with a current increasing 
consistently follows next stages: a) slow increasing at 
< j < j ca in the TAFF regime, where v a 1, b) sharp 
step-like increasing with a width of the order of j ca which 
corresponds to nonlinear transition from the TAFF to the 
FF regime, c) second stage of slow increasing for j > 2j ca 
which corresponds to the FF regime (see also item C of 
Sec. II). It follows from the expression for j ca (a,q) that 
an increasing of a and (or) q leads to a broadening of the 
step of the order of j ca and its shift to the larger current 
densities j w j oa . 

The anisotropy of fi(a) (see Eq. (13"8"|) and Fig. 4) can 
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FIG. 5: The current-angle dependence of the isotropic prob- 
ability function Ui(j,a). 



be divided into two types: simple ("external") which de- 
pends on cos 2 a, and more complex ("internal"), given 
by v a (a). The first (external) anisotropy stems from the 
"tensorial" a-dependence which exists also in the linear 
(TAFF and FF) regimes of the flux motion. The second 
(internal) is through the a-dependence of v a , which in 
the region of transition from the TAFF to the FF regime 
is substantially nonlinear (Eq. (j3"i2|) and Fig. 3). The ap- 
pearance of nonzero sin 2 a term in /, for a physically 
describes the guiding of vortices along the channels of the 
PPP in the presence of i-pins for the current densities 
3 ?5 jd{ct)- The influence of ;/ a -anisotropy on is dif- 
ferent for different values of the angle a (see Fig. 5). For 
a > 7r/4 the anisotropy of v a (a) does not influence the 
value of fi(a) because (f a /tana) 2 <C 1 in the expression 
f° r jd{ot)- On the contrary, for a <C 7r/4 the influence 
of a-pins on Vi (a) is most effective for that range of cur- 
rent density, where v\ > tan 2 a, due to the inequality 
tan 2 a <C 1. Thus, the i/j and v a as functions of the 
angle a at j = const behave themselves oppositely (see 
Figs. 3, 5): v± increases mo notonically with a-increasing, 
whereas v a - monotonically decreases. For j > j ca {a) 
and at small angles which meet the condition tan 2 a <C 1, 
the behavior of the V{ and v a qualitatively similar in a 
and opposite in q. 

In case where tan 2 a > 1, the i>i and v a behav- 
ior is qualitatively different and stems from the (a, in- 
dependences of the corresponding crossover current den- 
sities. In contradistinction to is a , the transition of vi 
from the TAFF to the FF depends weakly from a and q; 
it moves to the lower current densities with g-increasing 
for a > 7r/4 and moves to the higher ones for a -C 7t/4. 
In general, the v a behavior is more anisotropic than Vi 
behavior. The Vi anisotropy appears only in the TAFF 
regime, whereas v a anisotropy exists as in the TAFF, as 
well in the FF regime. And this anisotropy is greater in 
the current density as the angle a is greater. The Vi and 
v a transition width at a = const is defined by Si and e a 




FIG. 6: The current-angle dependence of the dimensionless 
even longitudinal anisotropic magnetoresistivity pf (j, a) for 
the value of the parameter e a = 1. 




FIG. 7: The current-angle dependence of the dimensionless 
even longitudinal magnetoresistivity pt(j, ol) for the value of 
the parameter e a = 1. 

parameters, respectively, and it increases for e% — > 1 and 



2. [j, a) -presentation of even magnetoresistivities. 

Now we are in a position to discuss the results of the 
presentation of Eq. (|25|) in the form of graphs. First 
we note that according to Eqs. (f!?5]) . the even resistive 
responses can be represented as the products of corre- 
sponding isotropic and anisotropic ^-functions. For this 
reason the graphical analysis of the pt(j> an d p\{j, ot), 
after the above-mentioned consideration of the vi (j, a) 
(see Fig. 5), can be reduced to the construction and anal- 
ysis of the pj a (j,a) and p\ a {j,a) graphs. 

Let us begin with a discussion of pf behavior (see 
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FIG. 8: The current-angle dependence of the dimensionless 
even longitudinal magnetoresistivity pt(j, a) for the value of 
the parameter e a = 0.01. 



Eq. (HU) and Fig. 6). For all a ^ due to the term 
sin 2 a in Eq. (118|) a critical current density j c exists only 
for direction, which is strictly perpendicular to the PPP 
(a = 0) (as it was shown in item E.l of Sec. II) and 
j c {a) = for any other direction (0 < a < it/2) due to 
the guiding of vortices along the channels of a washboard 
potential (see also Fig. 8 in Ref. 11). In the FF-regime the 
isotropization of the pjJ" Q arises due to the vortex slipping 
over the PPP channels. Thus at small angles a the v a 
function strongly influences the pjf, whereas for a — > 

7r/ 2 this influence is not so effective due to the external 
anisotropy, which is proportional to the sin 2 a term. 

Returning now to the consideration of the pt (j, a) 
graph we refer to the Eq. (|33|) . It is necessary to pay 
special attention to the TAFF behavior of these curves 
at small currents and temperatures, which follows from 
the full pinning of vortices by point-like pins. This be- 
havior is completely different (for a ^ 0) from the non- 
TAFF behavior of the corresponding graphs for the case 
of purely anisotropic pinning (see Fig. 8 in Ref. 11), 
which is provocated by the guiding of vortices along the 
channels of the PPP. At high current densities and (or) 
temperatures appears the FF regime, because the vortex 
motion transverse to the a-pins becomes substantial and 
longitudinal resistivity practically becomes isotropic. In 
these limiting cases the pt{j) magnitudes are equal to 
unity (Fig. 7). 

For the angles < a < it/2 the P^(j) behavior fol- 
lows substantially the properties of one multiplier. The 
qualitative behavior of these multipliers, depending on 
the j and a magnitude is very different as determined by 
different behavior of their crossover current densities jd 
and j ca - The priority of a sharp rise of the appearance 
U{ or v a functions depends on the competition between 
the crossover current densities jd an d j C ai respectively. 
That is why it may appear a "step" on some of the pt(j) 




FIG. 9: The current-angle dependence of the dimensionless 
even transverse anisotropic magnetoresistivity p± a (j,ot) 



curves (for q > 1 and a ^ 0,7r/2) when the next se- 
quence of the vortex motion regimes is realized: a) full 
i-pinning in the TAFF regime (0 < j < jd)l b) non- 
linear transition from the TAFF to the FF regime for 
z-pins (j > 2j ci ), c) practically linear the FF regime as 
a consequence of the guiding of vortices along the chan- 
nels of the washboard PPP (on the pjj" (j, a) surface one 

can see the horizontal sections at j w j ca , see Figs. 7, 8); 
d) nonlinear transition to the FF regime of vortex motion 
transverse to the a-pins for j > j ca and, at last, e) a free 
FF motion for j ^> j ca ■ 

With decreasing of the q the a)-e) corresponding re- 
gions along the current density axis j can overlap each 
other and a common nonlinear transition appears instead 
of b)-d) regions. For the limiting cases a — 0, tt/2, a 
guiding of vortices is absent and the pt(j) LT-behavior 
is simply related to the Vi and v a behavior. If parameter 
e a is decreasing, then the width of the transition of v a 
from the TAFF to the FF is also decreasing. Such en- 
hancement of the v a steepness leads to appearance of the 
minimum in a for the P|j~(i, a) graph (see Fig. 8). 

Now we pass to a discussion of the p\ a (j, a) and 
PjlUi a ) graphs. As it follows from Eq. (fT8]) . the p\ a < 
and has a minimum in a for all j = const. The p~\_ 
reaches its maximal magnitude for a ~ 7r/4 due to the 
factor sin a cos a and realization of guiding in the TAFF 
regime for j < j ca (see. Fig. 9a in Ref. 11 and Fig. 9). 
Therefore, the most favorable angle for its observation 
is near a = tt/A. In considered case the origin of this 
minimum has the same reason as a low (j, a)-behavior 
of the py (i) oi) curves in Fig. 7, namely it stems from 
existence of the TAFF regime for the point-like pins at 
small j-values. As is seen in Fig. 9, the position and 
the magnitude of this p^ a -minimum strongly depends on 
the a- value. It is very much pronounced for q > 1 and 
strongly suppressed for q < 1 by influence of the z-pins. 
With increasing of the current density j > j ca a position 
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FIG. 10: The current-angle dependence of the dimensionless 
even transverse magnetoresistivity p+ ( j, a) . Pay attention to 
the inverted direction of the axes in comparison with Fig. 9. 



of the minimum in a is shifting due to the competition 
of two multipliers in the p± a {j, ex.) expression Eq. (|18ll ): 
sin 2a is decreasing for a — > n/2, whereas (1 — v a (j, a)) is 
increasing with a-increasing for j — const, and decreas- 
ing for j-increasing for a — const due to the transition to 
the FF regime. For all a and current densities j > 2j ca 
the u a « 1, and for this reason p~^ a — > 0. The g-influence 
is defined by j C a(oi) and determines the region of appear- 
ance of a small value of the p\ a for the current densities 

J ~ 3ca- 

Since the pj_, according to Eq. is the product 

of the p\ a and i>i(fi), so this graph (see Fig. 10) can 
be reduced to the product of the graphs in Fig. 5 and 
Fig. 9. The transition from the TAFF to the FF regime is 
highly anisotropic in a; this causes a shift of the maximal 
p\_ (j, a) magnitude in the direction of a small angle a <C 
7r/4 for the j = const. That is why in view of i-pinning 
presence the p~^_(j,a), as distinct from p\Jj, ct), has the 
minimum both in a and in j. This statement follows from 
the fact that influence of i-pinning leads to the p\_ — > 
for < j < jci(d) due to the i/jCl, For the current 
densities j > j ca (a) the p~\_ (j, a) behavior is determined 
exclusively by the above-mentioned p\ a {j,a) behavior. 



3. {j, a) -presentation of odd magnetoresistivities. 

Before following discussion of the odd resistive re- 
sponses we should remind the reader about the bump-like 
behavior of the current and temperature dependence of 
the v~ functions (see Figs. 6 and 7 in Ref. 11 ), because 
Z/r and v~ functions, as it follows from Eqs. (f55|) - ([5rj)) . 
give an important contribution to the odd responses. The 
v~(j) and v~~(0) curves for the case of e < 1 in fact 
are proportional to the derivatives of the corresponding 
v + {j) and v + {0) curves, which have a step-like behav- 




FIG. 11: The current-angle dependence of the function 




FIG. 12: The current-angle dependence of the function 



ior as a function of their arguments (see Ref. 11 for the 
detailed discussion of this point and Eq. (fTT|) in this pa- 
per). As the p7 and p^ resistivities given by Eqs. (|26|) - 
(|2"T)) arise by virtue of the Hall effect, their characteristic 
scale is proportional to e <C 1, as for Eqs. (|19p for purely 
anisotropic pins. 

The position of the characteristic peak in the v~ and 
v~ functions is different for q ^ 1, because parameter q 
determines the anisotropy of the critical current densities 
for i- and a- pins. So, if q is not very close to the unity, 
the position of the i- and a- peaks cannot coincide, and 
in this case the current and temperature odd resistive 
dependences pj ± can have a bimodal behavior. For the 

pjj" curves such dependences will correspond to existence 
of the resistive "steps" on the pjj~ curves (see Fig. 7), 
because for e <C 1 we can consider the pjj" dependences 
as derivatives of the pi curves. From this viewpoint it is 
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FIG. 13: The current-angle dependence of the odd longitudi- 
nal magnetoresistivity p7(j,a). 

easy to understand the previous assertion in item E.2 of 
Sec. II that pjj" includes two terms (every proportional 

to the v~ and v~ , respectively) with similar signs. 

Now we will discuss the v~ and za~ as a function of 
(J, a) and the parameter q in detail. Really, due to the 
smallness of the Hall constant, the v~ and v~ tend to zero 
in the regions of the linear TAFF and FF regimes of the 
v a and Vi function, respectively. The v~ and v~ functions 
have a sharp peak (see Fig. 11, 12) in the region of sharp 
change of the v a and Vi increasing (for j w j ca or j « 
j c i, respectively). With a- and q-increasing the width 
and the height of the v~ maximum also increases with 
simultaneous shift of the maximum to the higher current 
densities due to the relation j ca (a) ~ q/cosev. The v~ 
peak is located in the angle range < a < ir/4, which 
corresponds to a change of the angular dependence of the 
crossover current density j C i(a) from the angles a > tt/4 
to the angles a <C 7r/4 (see. Fig. 12). The v~ maximum 
shifts to a smaller current densities with q-increasing due 
to the jd(oi) w l/vaqcosa,. The magnitudes of the v~ 
and v~ are compete by an order of magnitude for < a < 
7r/4 and all q- values which satisfy a condition j ca « j c ;. 

Now let us discuss a graphical presentation of the 
Eq. (f3"5|) , which can be represented as pjj~ = B\ + £> 2 , 

where B\ = ^ pt a , and B2 — PfV~Vi cos 2 a. Taking into 
account that every factor in the B\ and B2 is positive 
(see Figs. 5, 6, 11, 12), we can conclude that p^ > for 
all values of the j, a, q. 

Proceeding to the analysis of the B\ and B2 (j, a, q)- 
behavior in details we consider first those limiting cases 
in which a- or i- pinning is dominant i.e. Vi « 1 or 
i/„ w 1, respectively. If a-pinning is dominant (i.e. for 
q 3> 1), then — > 0, and Eg. ([33]) has the form p^ w 

/OjT = i/~ cos 2 a. For the opposite case (i.e. for q < 1), 

conversely, v~ -> 0, and w ^ r P|f Q - The Pj^O':") 
graph presentation is especially simple because it may 




FIG. 14: The current-angle dependence of the odd trans- 
verse magnetoresistivity p~J_(j, a). The characteristic mini- 
mum (which is shown by the arrow) is in the region < a < 
7r/4 and j ~ j c i(a). The minimum is shown as two neighbor- 
ing minimums due to the step-like behavior of the calculation. 
Pay attention to the inverted direction of the axes in compar- 
ison with Fig. 13. 



be depicted with the aid of Figs. 5, 11, 12. 

In the general case, i.e. for q « 1, we should con- 
sider the B\ and B2 separately because dominant type 
of pinning is absent. The B\ is proportional both u~ , 
which is nonzero for < a < tt/4 and j{a) ks j c i{oL,q) 
(see Fig. 12), and the factor /?jj~ o with a graph, shown in 
Fig. 6. As a result, the B\ has a sharp maximum for the 
< a < 7r/4 and j(a) w jd(oi, <?)■ The second term B2 is 
proportional both the factor v~ and the factor z/j cos 2 a. 
The contribution of the first factor is maximal for a sa 
7r/2 and current densities j(a) ~ j ca (a,q), whereas the 
Vi cos 2 a contribution is maximal for < a < 7r/4 and 
j(oi) j C i{oi, q)- Therefore these factors compete so that 
the resulting maximum of B2 shifts from a « n/2 to the 
tt/4: < a < tt/2. It is relevant to note that the con- 
dition j c i (a) < j C a{&) for q « 1 and 7r/4 < a < tt/2 
is always fulfilled. That is why the maximal contribu- 
tion of the B2 is realized for j(a) « j ca (a.,q) because 
in this region of the current densities the v% — * 1 for 
7r/4 < a < tt/2. Therefore the pjj" behavior is determined 
mainly by the B2 behavior, and the B\ contribution is 
essential for < a < 7r/4 and j w jci (a). 

The /Oj(j, a) dependence is the most complicated. For 
the sake of simplicity the analysis we represent the p^_ as 
asumpj = pf[Ai+(A 2 +A 3 )sm2a], where A 1 = nev a ff, 
A 2 = v~Vi/2, A 3 = —v^{l — v a )/2. First we consider 
the limiting cases of purely isotropic or anisotropic pin- 
ning [y a — » 1 or V; — > 1, respectively). For i-pinning 
we have p^_ = pfneuf, from which follows (Vinokur et 
al. 17 ) a scaling relation p± ~ p 2 . For the case of purely 

anisotropic pinning = Pf{ nev a + (v~ sin2a)/2}, and 
the scaling relation is p± ~ pii (see also Ref. 16). 
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Now we consider every term in the p J (j, a) in detail. 
The Ai contribution can be reduced in fact to the mul- 
tiplication of the graph in Fig. 3 by the graph in Fig. 5 
squared; the result is essentially nonzero for j > j ca (&, q)- 
The A 2 contribution was described above (see the B 2 
term in the p^ without taking into account the cos 2 a 
anisotropy). Note also that both terms (A\ and A 2 ) are 
positive for ne > 0. The A3 behavior is of great interest 
because the A 3 < for ne > 0. Let us consider the cases 
q > 1 and q < 1, which correspond to the a-, or i-pinning 
domination, respectively. Then, for a < tt/4: 

a) for q < 1 we have jd(a) > j ca (a) and the sharp 
maximum of the v~ is suppressed by the factor (1 — u a ) — > 
0. As a result, the A 3 contribution can be ignored. 

b) for q > 1 the opposite inequality follows, i.e. 
jd(a) < jca(a)- Then for j w j ci (a) the A 3 term is 
dominant because f a < 1 and i/T" — > ne in this (j, a)- 
rcgion (see Fig. 3 and Fig. 5). As a result, the p±(j, a, q) 
change the sign for j w j C i(a) and < a < 7r/4. Since 
the scale of the <C fj, the amplitude of the minimum 
is small in comparison with the pj magnitude. 

Thus, a competition of the a- and i-pinning leads 
to the qualitatively important conclusion that the p~]_ 
can change its sign at a certain range of (a,j, g)-values, 
namely for j w j ci (a, j, g), < a < tt/4, and g > 1. 

D. Resistive response in a rotating current scheme. 

1. Polar diagram. 

An experimental study of the vortex dynamics in 
YBa2Cu307_a crystals with unidirectional twin planes 
was recently done using a modified rotating current 
scheme 4,5 . In that scheme it was possible to pass current 
in an arbitrary direction in the ab plane of the sample 
by means of four pairs of contacts placed in the plane 
of the sample. Two pairs of contacts were placed as in 
the conventional four-contact scheme, and the other two 
pairs were rotated by 90° with respect to the first (see 
the illustration in Fig. 1 of Ref. 4). By using two cur- 
rent sources connected to outer pair of contacts, one can 
continuously vary the direction of the current transport 
in the sample. By simultaneously measuring the voltage 
in the two directions, one can determine directly the di- 
rection and magnitude of the average velocity vector of 
the vortices in the sample as a function of the direction 
and magnitude of the transport current density vector. 
This made it possible to obtain the angular dependence 
of the resistive response on the direction of the current 
with respect to the pinning planes on the same sam- 
ple. The experimental data 4 ' 5 attest to the anisotropy 
of the vortex dynamics in a certain temperature inter- 
val which depends on the value of the magnetic field. A 
rotating current scheme was used to measure the po- 
lar diagrams of the total magnetoresistivity p(a), where 
p = (p 2 , + Py) 1 ^ 2 is the absolute value of the magnetore- 
sistivity, p x and p y are the x and y components of the 




FIG. 15: Series of graphs of the function p(a) for a sequence 
of the parameter j: 0.63 (1), 0.65 (2), 0.75 (3), 1.00 (4), 1.50 
(5), 1.92 (6), 2.00 (7), 2.50 (8), 4.00 (9), 20.0 (10) for e a = 1. 



magnetoresistivity in an xy coordinate system, and a is 
the angle between the current direction and the oy axis 
(parallel to the channels of the a-pinning centers). In the 
case of a linear anisotropic response the polar diagram of 
the resistivity is an ellipse, as can easily be explained. In 
the case of a nonlinear resistive response the polar dia- 
gram of the resistivity is no longer an ellipse and has no 
simple interpretation. 

In this subsection we carry out a theoretical analysis 
of the polar diagrams of the magnetoresistivity p in the 
general nonlinear case in the framework of a stochastic 
model of a + i pinning. This type of angular dependence 
p(a) is informative and convenient for theoretical analy- 
sis. For a sample with specific internal characteristics of 
the pinning (such as q, e a , £j, and k) at a given temper- 
ature and current density the function p(a) is contained 
by the resistive response of the system in entire region of 
angles a and makes it possible to compare the resistive 
response for any direction of the current with respect to 
the direction of the planar pinning centers. In addition, 
in view of the symmetric character of the p(a) curves, 
their measurements makes it possible to establish the 
spatial orientation of the system of the planar pinning 
centers with respect to the boundaries of the sample if 
this information is not known beforehand. 

Now for analysis of the p(a) curves we imagine that 
vector j is rotated continuously from an angle a = it/ 2 
to a = 0. The characteristic form of the p(a) curves will 
obviously be determined by the sequence of dynamical 
regimes through which the vortex system passes as the 
current density vector is rotated. By virtue of the sym- 
metry of the problem, the p(a) curves can be obtained 
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FIG. 16: Series of graphs of the function p(a) for a sequence 
of the parameter j: 0.63 (1), 0.65 (2), 0.75 (3), 1.00 (4), 1.50 
(5), 1.92 (6), 2.00 (7), 2.50 (8), 4.00 (9), 20.0 (10) for e a = 0.1. 



in all regions of angles a from the parts in the first quad- 
rant. 

We recall that in respect to the two systems of pin- 
ning centers it is possible to have the linear TAFF and 
FF regimes of vortex dynamics and regimes of nonlinear 
transition between them. The regions of nonlinear tran- 
sitions are determined by the corresponding values of the 
crossover current densities j C i(a, q) and j ca (&, <?)• 

Now let us consider the typical p(a) dependences which 
are presented in Fig. 15 and 16 for a sequence of a current 
density magnitude. We remind that the polar diagram 
graphs represented below are constructed, as the previous 
graphs in Figs. 3-7, 9-14 for the next values of param- 
eters: q — 1.6 (i.e. for the case with dominant a-pins), 
k = 0.25, 6 = 0.003, e = 0.01, = 0.1, e a = 1 (Fig. 15), 
and e a = 0.1 (Fig. 16). Note that p(a) is the product 
of two multipliers: one is the Vi(fi) dependence, which 
was earlier studied in Fig. 4 of item C.l of Sec. Ill, and 

other is the \J sin 2 a + cos 2 a factor, which qualitative 
behavior is close to the pt a (j, (x) dependence (see Fig. 6 
in item C.2 of Sec. III). 

Let us analyze the p{a) behavior for the series of val- 
ues of the current density j. When the angle a changes 
from to tt/2 the function p(a) grows monotonically 
from p(0) = v a {j I 'q)vi{jqv a {j / 'q)) to p(n/2) = v l (]q). 
In Fig. 15 curves 1-6 of the function p(a) have the shape 
of the 8-figure drawn along ox-axis (strongly elongated 
for the curves 1,2). 

This anisotropy can be determined by the relation of 
the magnitudes of the half-axis at the direction a = n/2 
to the transverse half-axis for any fixed magnitude of the 



current density. The curves 1-6 of the p = p{a) graph 
has the 8-form elongated along ox-axis. It is caused by 
the step-like behavior of the function, corresponding 
for the curves 1, 2 to the crossover from the TAFF to the 
FF regime. That is why the magnitude of the p(ir/2) for 
the curve 2 is rather greater than for the first one. With 
a-increasing the ^-function is in the TAFF-region (see 
Fig. 5), which provocates the p(a) <C 1 in the case where 
the condition j < j c i{dt) is satisfied. Therefore, with j- 
increasing the magnitude of the angle a, which separates 
the TAFF and the FF regions of the ^-function at a fixed 
value of the current density, decreases to the a<l, 

As the Vi(j,a = tt/2) is in the FF region (i.e. 
3 ^ jci(ct = tt/2)), so the anisotropy of the 8-curvc 
decreases for curves 3-6. The p = p(a) behavior of the 
curves 5-6 is more isotropic in the region a <C 7r/4 than 
behavior of the curves 1-4. If the condition j > j ca (0) is 
satisfied, an appearance of the nonzero resistance in cor- 
responding region follows. Its magnitude is smaller than 
p(tt/2) for the curves 7, 8, 9 and practically is equal to 
the p(n/2) for the curve 10. Note, that for the a = 0,7r 
and j ca < j < 3j ca /2 one can see the minimum, which 
decreases with j-increasing and disappears in the case 
where the condition j > 3j ca /2 is satisfied. So, for large 
magnitudes of the current densities the p(a) behavior be- 
comes more isotropic. 

It is necessary to pay attention for the p = p{a) behav- 
ior in the case where s a = 0, 1 (see Fig. 16) for the same 
series of the magnitudes of the current densities. The 
behavior of the curves 6, 7, 8, 9 differs from the above- 
mentioned case, but the behavior of the curves 1-5, 10 
retains the same. This fact is caused by the influence of 
the parameter e a on the Vi behavior only in the area of its 
sharp step-like behavior at the j ~ j ca (a). Note, that the 
v-i contribution is dominant in the region a < 1 as well 
as the above-mentioned anisotropy of the p\ a (j, a) (sec 
Fig. 5 and Fig. 7). As decreasing of the e a causes the 
more narrow crossover from the TAFF- to the FF-regimc, 
the Vi{a) has a minimum at fixed magnitude of the cur- 
rent density. The magnitude of this minimum decreases 
with the j-increasing and the minimum shifts from the 
a ~ 7r/4 to the a ~ 7r/2. The influence of the parameter 
q acts on the crossover current densities jd and j ca only 
quantitatively, but does not change an evolution of the 
curves 1-10 qualitatively. 



2. Qe (a) -dependence. 

Let us examine theoretically in our model a new type 
of the experimental dependence, recently studied 4 for 
6_E(a), where 9^ is the angle between j-vector and 
the electric field vector E measured at fixed values of 
the current density and temperature. Taking into ac- 
count that in the xy coordinate system the magnetore- 
sistivity components are p x — p xx since = Vi(Fi) sina, 
p y = p yy cos (x = v 'j (.Fj ) v ' a (f a) cos a, we obtain the follow- 
ing simple relation: tan8£(a) = p x /Py = tana/ 'is a (F a ), 
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FIG. 17: Series of graphs of the function O e (ct) for a sequence 
of the parameter j: 1 (1), 1.7 (2), 2.2 (3), 3.5 (4), 20 (5) for 
T = SK. 
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FIG. 18: Series of graphs of the function j c (ct) for the pa- 
rameter pairs: E c = 0.002, 5=1.6 (1); E c = 0.002, g=3 (2); 
Ec = 2, 9 =1.6 (3). 



Oe (a) — arctan(tana/^ a (i 7 ' a )). 



(39) 



Note, that the vi term, describing the i-pinning, is ab- 
sent in Eq . (f3"9"|) . Then it follows from the latter that 
the v a (j, a , 0) — tana/ tan Qe{j, ol,9), i.e. the v a {j, a, 9) 
function can be found from the experimental dependence 
Qe{c(). Unfortunately, the dependence @e((x) for the se- 
ries of the temperature values was experimentally found 4 
so far only for the FF-regime (see. Fig. 2 in Ref. 4). The 
&eU, ct) dependence is presented in Fig. 17. It shows 
all changes in the Qe(j, ot) behavior also for the TAFF- 
regime. 

Let us analyze the Eq. (|39[) in detail. The O e (j, ct) is 
the odd function of the angle a, and its magnitude in- 
creases monotonically with the a-increasing for all values 
of the j due to the monotonical decreasing of the v a (j, a) 
function (see. Fig. 17). It follows from Eq. (|39|) that the 
period of the function ©£j(a) is equal to it. One more 
important limiting case is realized for v a w 1, which cor- 
responds to the limit of isotropic pinning. Depending on 
the inequality between the j magnitude and the crossover 
current density j ca {oi) ~ qj cosev, one can separate two 
regions where the @e(J, ol) behavior is qualitatively dif- 
ferent. If A is the argument of the arctangent function in 
Eq. (|39p. then in that region j, a, q, where the inequality 
j <L jca{ot) is true (the FF regime for u a (J, a), see also 
Fig. 3), the magnitude of the 9g ~ i as i C 1. And 
for the case j < j ca (the TAFF regime of the v a (j,a)) 
the value Q E ~ n/2 - A^ 1 , as A > 1. 

Note, that the parameter e a influences the 6 e {j: ce) by 
changing the character of the step-like crossover of the 
v a (J,a) (the smaller the e a , the sharper the crossover). 



The value of the parameter q, as well as above-mentioned, 
determines the magnitude of the j C a( a ) (and, therefore 
the position of the boundaries in j of the regions of quite 
different 0£(a) behavior) at fixed a. 



3. Critical current density anisotropy. 

Under the critical current density we mean the current 
density, which corresponds to the electric field strength 
on the sample E = 1/iV '/cm. Let us determine the j c (ct) 
behavior graphically by crossing the = jp^(j) graph 
and the plain E = E c in the polar coordinates. For all 
angles a the point of crossing for these graphs determines 
the critical current density magnitude for the defined di- 
rection, and the crossing line of the graphs presents the 
dependence j c (&)- 

Let us remind the reader that as in above-mentioned 
sections, in the nonlinear law E^ = jp^(j) we measure 
j and p in the values of the jo = cUo/^odh and Pf — 
p n B/B C 2, respectively. That is why the E magnitude we 
have to measure in the Eq = joPf. As well as in item 
B of Sec. Ill we use the data from Ref. 8, where for the 
niobium samples p n w 5, 5 ■ 10 -6 Ohm-cm, B w 150 Gs, 
B c2 w 17 kGs, p f w 5 ■ lCT 8 Ohm-cm, U = 2500K, and 
d=2.5- 1(T 6 cto. 

Therefore, E w 6 • 10 _4 U/cm, and for E c = lpV/cm 
we have to cross the dimensionless pt(j) '3 graph by the 
plain E « 0.002. 

Now we will discuss the j c {ot) as a function of a, q, E Cl 
and e in detail. The j' c (cv)-anisotropy can be determined 
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by the relation of the magnitudes of the half-axis at the 
direction a = to the transverse half-axis for any fixed 
magnitude of the parameters q, E c . The j c (o) decreases 
monotonically from j c (0) with a-increasing and has a 
minimum for a = n/2. It is caused by the fact that, 
as it was shown in item C.l of Sec. II, the a-pinning 
(with high values of the j ca for q > 1) does not influence 
the i-pinning for a = tt/2. Therefore, the inequality for 
the crossover current densities jd(a) < j ca {ct) for q > 1 
leads to the corresponding inequality for the critical cur- 
rent densities j c (0) < j c (7r/2). 

The q influences the j c {a) behavior (as in item D.l 
of Sec. II) only quantitatively: with g-increasing the 
ratio j c (0)/j c (7r/2) grows and visa versa. It is caused by 
the j c (0)-increasing and j c (7r/2)-decreasing due to the a- 
behavior of the corresponding crossover current densities 
jci(a) and jca(ck). The smaller the e , the sharper the 
crossover between the j c (a) regions of slowly and quickly 
decreasing as a function of the a. With E'c-increasing 
the nonlinear law = p^(j)j is satisfied for the larger 
values of the current density. 

That is why with a-increasing from to a* values 
(for which the condition tan 2 a* <C v^(j,a*) is satis- 
fied) the f a -function is in the FF regime and j c (ct) de- 
creases slowly. When the condition a > a* is true the 
i^-function has a step-like crossover from the FF to the 
TAFF regime and j c (a) decreases quickly. 

So, the a* behavior as a function of the parameters 
q and E c is qualitatively different: it increases with E c - 
increasing and decreases with (/-increasing. On the in- 
crease of the E c by the several orders of magnitude the 
j c (a) curve degenerates into a circumference due to the 
isotropization of the j C i(a) and j C a{oL) behavior for the 
high j-values. Otherwise, with -E^-decreasing the j c {o) 
curve degenerates into a narrow loop, because the j C i{pt) 
and j ca (a) behavior for a small j is very anisotropic. 



IV. CONCLUSION. 

In the present work we have theoretically examined the 
strongly nonlinear anisotropic two-dimensional single- 
vortex dynamics of a superconductor with coexistence 
of the anisotropic washboard PPP and isotropic pinning 
potential as function of the transport current density j 
and the angle a between the direction of the current and 
PPP planes at a fixed temperature 9. 

The experimental realization of the model studied here 
can be based on both naturally occurring 2-5 and artifi- 
cially created 6-8 systems with i + a pinning structures. 
The proposed model has made it possible for the first 
time (as far as we know) to give a consistent description 
of the nonlinear anisotropic current- and temperature- 
induced depinning of vortices for an arbitrary direction 
relative to the anisotropy of the washboard PPP. In the 
framework of this model one can successfully analyze the- 
oretically certain observed resistive responses which are 
used for studying anisotropic pinning in a number of new 



experimental techniques 4 (the polar diagram of p(a), the 
9e(a) curve described by formula Eq. (|3H)) ) as well as 
new Hall responses specific for the i + a pinning problem. 

A quantitative description of the anisotropic nonlin- 
ear resistive properties of the problem under study is 
done in the framework of the stochastic model on the 
basis of the Fokker-Planck approach. The main nonlin- 
ear components of the problem are the anisotropic v a (F a ) 
and isotropic Vi(Fi) probability functions for the vortices 
to overcome the potential barriers of a- and i-pinning 
centers under the action of anisotropic motive forces F a 
and Fi, respectively. The latter include both the "ex- 
ternal" parameters j,a,9 and the "internal" parameters 
q,Ei,e a which describe the intensity and anisotropy of 
the pinning. As can be seen from Eqs. (f3"3"| - (j3l))) , the 
magnetoresistivities Py ot, 0) are, in general, nonlin- 
ear combinations of the experimentally measured vi and 
v a functions (vi can be measured independently from the 
= tt/2), see Eq. (|33|) and v a - from the 9_e(o:), 
see Eq. (15^)) ). 

Therefore, the nonlinear (in j) resistive behavior of 
the vortex system can be caused by factors of both an 
anisotropic and isotropic pinning origin. It is important 
to underline that whereas the structure of the v a {F a ) and 
F a is the same as for purely a-pinning problem, the struc- 
ture of the Vi(Fi) and Fi is strongly different from the 
structure of the purely i-pinning problem due to the fact 
that Fi , as motive force of the (i + et)-problem, is nonlin- 
ear and anisotropic (see Eqs. (|3"7) - (|3"5)0 and Figs. 3, 4, 5). 

Two main new features appear due to the introduc- 
tion of the isotropic i-pins into the initially anisotropic 
a-pinning problem. First, unlike the stochastic model 
of uniaxial anisotropic pinning studied previously 10,11 , 
where the critical current density j c is indeed equal to 
zero for all directions (excepting a = 0) due to the guid- 
ing of vortices, in the given i + a model the anisotropic 
critical current density j c (a) exists for all directions be- 
cause i-pins "quench" the guiding of vortices in the 
limit (j, 6) — > 0. Second, the Hall resistivity response 
functions p~]_(j,a) can have a change of sign in a cer- 
tain range of (J, a, q) (at fixed dimensionlcss Hall con- 
stant e = olh/t} and the dimensional Hall conductivity 
5 = ne/pf), whereas the sign of the p~^(j,a) does not 
change. 

It should be noted that recently 8 the nonlinear (in 9) 
anisotropic longitudinal and transverse resistances of Nb 
films deposited on facetted sapphire substrates were mea- 
sured at different angles a between j and facet ridges in a 
broad range of temperature and relatively small magnetic 
field H. The experimental data were in good agreement 
with the theoretical model described here. The measured 
pt(9,a) dependences can be fitted using the probabil- 
ity functions v a and in the form proposed here (see 
Eq. (|32p ) with the anisotropic and isotropic pinning po- 
tential given by Eq. (|3"Tj) . The periods and depths of 
the potential wells were estimated from the experimental 
data 8 and were used here (see Sec. Ill) for the theoret- 
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ical analysis of different types of nonlinear anisotropic 
(J, a)-dependent magnetoresistivity responses, given by 
Eqs. ([33 ]) - ([56j) . in the form of graphs (see Figs. 3-18). 
Whether these theoretical results can explain a new por- 



tion of the (J, a)-dependent i+a resistivity data measured 
at 9 = const (in particular, for the samples investigated 
earlier 8 at small current densities) remains to be seen. 
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